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Abstract 

h— •} In previous works, we showed that the internal DLA cluster on 1 d with t 

particles is almost surely spherical up to a maximal error of 0(log t) if d = 2 and 
ff) 0(v / logi) if d > 3. This paper addresses "average error": in a certain sense, 

the average deviation of internal DLA from its mean shape is of constant order 
'-^j ' when d = 2 and of order r 1_d / 2 (for a radius r cluster) in general. Appropriately 

normalized, the fluctuations (taken over time and space) scale to a variant of 

the Gaussian free field. 
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1 Introduction 

1.1 Overview 

We study scaling limits of internal diffusion limited aggregation ("internal DLA"), 
a growth model introduced in [MD86, DF91]. In internal DLA, one inductively 
constructs an occupied set At C Z rf for each time t > as follows: begin with 
Aq = and A\ = {0}, and let A t +i be the union of A t and the first place a random 
walk from the origin hits "L d \ Af. 

The purpose of this paper is to study the growing family of sets At- Following 
the pioneering work of [LBG92], it is by now well known that, for large t, the 
set At approximates an origin-centered Euclidean lattice ball B r := B r (0) n 7j d 
(where r = r(t) is such that B r (0) has volume t). The authors recently showed 
that this is true in a fairly strong sense [JLS09, JLSlOa, JLSlOb]: the maximal 
distance from a point where l^ t — lB r is non-zero to dB r (0) is a.s. O(logt) if d = 2 
and 0(ylog f) if d > 3. In fact, if C is large enough, the probability that this 
maximal distance exceeds Clogt (or CyTogt when d > 3) decays faster than any 
fixed (negative) power of t. Some of these results are obtained by different methods 
in [AG 10a, AG 10b]. 

This paper will ask what happens if, instead of considering the maximal distance 
from dB r (0) at time t, we consider the "average error" at time t (allowing inner and 
outer errors to cancel each other out). It turns out that in a distributional "average 
fluctuation" sense, the set At deviates from B r (0) by only a constant number of 
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lattice spaces when d = 2 and by an even smaller amount when d > 3. Appropriately 
normalized, the fluctuations of At, taken over time and space, define a distribution 
on M. d that converges in law to a variant of the Gaussian free field (GFF): a random 
distribution on M. d that we will call the augmented Gaussian free field. (It can 
be constructed by defining the GFF in spherical coordinates and replacing variances 
associated to spherical harmonics of degree t by variances associated to spherical 
harmonics of degree £ + 1; see §1.5.) The "augmentation" appears to be related to a 
damping effect produced by the mean curvature of the sphere (as discussed below). 1 
To our knowledge, no central limit theorem of this kind has been previously 
conjectured in either the physics or the mathematics literature. The appearance 
of the GFF and its "augmented" variants is a particular surprise. (It implies that 
internal DLA fluctuations — although very small — have long-range correlations 
and that, up to the curvature-related augmentation, the fluctuations in the direction 
transverse to the boundary of the cluster are of a similar nature to those in the 
tangential directions.) Nonetheless, the heuristic idea is easy to explain. Before we 
state the central limit theorems precisely (§1.3 and §1.4), let us explain the intuition 
behind them. 

Write a point x G M. d in polar coordinates as r9 for r > and 9 on the unit sphere. 
Suppose that at each time t the boundary of A t is approximately parameterized by 
rt{9)9 for a function r% defined on the unit sphere. Write 

r t {e) = {t/u d ) l ' d + Pt {e) 

where u;^ is the volume of the unit ball in M. d . The pt(9) term measures the deviation 
from circularity of the cluster At in the direction 9. How do we expect pt to evolve 
in time? To a first approximation, the angle at which a random walk exits At is a 
uniform point on the unit sphere. If we run many such random walks, we obtain 
a sort of Poisson point process on the sphere, which has a scaling limit given by 
space-time white noise on the sphere. However, there is a smoothing effect (familiar 
to those who have studied the continuum analog of internal DLA: the famous Hele- 
Shaw model for fluid insertion, see the reference text [GV06]) coming from the fact 
that places where pt is small are more likely to be hit by the random walks, hence 
more likely to grow in time. There is also secondary damping effect coming from 
the mean curvature of the sphere, which implies that even if (after a certain time) 
particles began to hit all angles with equal probability, the magnitude of pt would 
shrink as t increased and the existing fluctuations were averaged over larger spheres. 

The white noise should correspond to adding independent Brownian noise terms 
to the spherical Fourier modes of pt- The rate of smoothing/damping in time should 
be approximately given by Apt for some linear operator A mapping the space of func- 
tions on the unit sphere to itself. Since the random walks approximate Brownian 

1 Consider continuous time internal DLA on the half cylinder (Z/mZ) d_1 x Z+, with particles 
started uniformly on (Z/mZ) d_:L x {0}. Though we do not prove this here, we expect the cluster 
boundaries to be approximately flat cross-sections of the cylinder, and we expect the fluctuations 
to scale to the ordinary GFF on the half cylinder as m — > oo. 
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motion (which is rotationally invariant), we would expect A to commute with or- 
thogonal rotations, and hence have spherical harmonics as eigenfunctions. With 
the right normalization and parameterization, it is therefore natural to expect the 
spherical Fourier modes of pt to evolve as independent Brownian motions subject 
to linear "restoration forces" (a.k.a. Ornstein-Uhlenbeck processes) where the mag- 
nitude of the restoration force depends on the degree of the corresponding spherical 
harmonic. It turns out that the restriction of the (ordinary or augmented) GFF on 
W 1 to a centered volume t sphere evolves in time t in a similar way. 

Of course, as stated above, the "spherical Fourier modes of pt" have not really 
been defined (since the boundary of A t is complicated and generally cannot be 
parameterized by rt{9)9). In the coming sections, we will define related quantities 
that (in some sense) encode these spherical Fourier modes and are easy to work 
with. These quantities are the martingales obtained by summing discrete harmonic 
polynomials over the cluster At. 

The heuristic just described provides intuitive interpretations of the results given 
below. Theorem 1.3, for instance, identifies the weak limit as t — > oo of the internal 
DLA fluctuations from circularity at a fixed time t: the limit is the two-dimensional 
augmented Gaussian free field restricted to the unit circle dB\(0), which can be 
interpreted in a distributional sense as the random Fourier series 
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where for k > and (3k for k > 1 are independent standard Gaussians. The 
ordinary two-dimensional GFF restricted to the unit circle is similar, except that 
y/k + 1 is replaced by \fk. 

The series (1) — unlike its counterpart for the one-dimensional Gaussian free 
field, which is a variant of Brownian bridge — is a.s. divergent, which is why we use 
the dual formulation explained in §1.4. The dual formulation of (1) amounts to a 
central limit theorem, saying that for each k > 1 the real and imaginary parts of 
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converge in law as r — > oo to normal random variables with variance 2 (k+i) ( an d 
that Mj and Mk are asymptotically uncorrelated for j ^ k). See [FL10, §6.2] for 
numerical data on the moments Mk in large simulations. 



1.2 FKG inequality statement and continuous time 

Before we set about formulating our central limit theorems precisely, we mention a 
previously overlooked fact. Suppose that we run internal DLA in continuous time 
by adding particles at Poisson random times instead of at integer times: this process 
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we will denote by A T {t) (or often just At) where T(t) is the counting function for 
a Poisson point process in the interval [0, t] (so T(i) is Poisson distributed with 
mean t). We then view the entire history of the IDLA growth process as a (random) 
function on [0, oo) x Z ', which takes the value 1 or on the pair (t, x) accordingly as 
x G Ax(t) or x ^ Aj<ny Write O for the set of functions / : [0, oo) x Z d — > {0, 1} such 
that f(t,x) < f(t',x) whenever t < t' , endowed with the coordinate-wise partial 
ordering. Let F be the distribution of {A T (t)}t>0i viewed as a probability measure 
on $7. 

Theorem 1.1. (FKG inequality) For any two increasing functions F,G G L 2 (Q,P), 
the random variables F({A T ^}t>o) and G({A T (t)}t>o) are nonnegatively correlated. 

One example of an increasing function is the total number jf=Arpu\C\X of occupied 
sites in a fixed subset X C Z d at a fixed time t. One example of a decreasing 
function is the smallest t for which all of the points in X are occupied. Intuitively, 
Theorem 1.1 means that on the event that one point is absorbed at a late time, 
it is conditionally more likely for all other points to be absorbed late. The FKG 
inequality is an important feature of the discrete and continuous Gaussian free fields 
[She07], so it is interesting (and reassuring) that it appears in internal DLA at the 
discrete level. 

Note that sampling a continuous time internal DLA cluster at time t is equivalent 
to first sampling a Poisson random variable T with expectation t and then sampling 
an ordinary internal DLA cluster of size T. (By the central limit theorem, \t — T\ 
has order y/t with high probability.) Although using continuous time amounts to 
only a modest time reparameterization (chosen independently of everything else) it 
is aesthetically natural. Our use of "white noise" in the heuristic of the previous 
section implicitly assumed continuous time. (Otherwise the total integral of pt would 
be deterministic, so the noise would have to be conditioned to have mean zero at 
each time.) 

1.3 Main results in dimension two 

For x G Z 2 write 

F(x) := mf{t : x G A T(t) } 

and 

L(x) := \/F(x)/n- \x\. 

In words, L(x) is the difference between the radius of the area t disk — at the time 
t that x was absorbed into At — and \x\. It is a measure of how much later or 
earlier x was absorbed into At than it would have been if the sets At{+) were exactly 
centered discs of area t. By the main result of [JLSlOa], almost surely 

L{x) 

nmsup : — - < oo. 

x&? log kl 
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(a) 



(b) 



Figure 1: (a) Continuous-time IDLA cluster Ax{t) for t = 10 5 . Early points (where 
L is negative) are colored red, and late points (where L is positive) are colored blue, 
(b) The same cluster, with the function L(x) represented by red-blue shading. 

The coloring in Figure 1(a) indicates the sign of the function L(x), while Fig- 
ure 1(b) illustrates the magnitude of L(x) by shading. Note that the use of contin- 
uous time means that the average of L(x) over x may differ substantially from 0. 
Indeed we see that — in contrast with the corresponding discrete-time figure of 
[JLSlOa] — there are noticeably fewer early points than late points in Figure 1(a), 
which corresponds to the fact that in this particular simulation T(t) was smaller 
than t for most values of t. Since for each fixed x £ 1? the quantity L(x) is a 
decreasing function of At(x), the FKG inequality holds for L as well. The positive 
correlation between values of L at nearby points is readily apparent from the figure. 

Identify IR 2 with C and let Hq be the linear span of the set of functions on 
C of the form Ke(az k )f(\z\) for a G C, k € Z>o, and / smooth and compactly 
supported on M>o- The space Ho is obviously dense in L 2 (C), and it turns out to 
be a convenient space of test functions. The augmented GFF (and its restriction to 
dBi(0)) will be defined precisely in §1.4 and §1.5. 

Theorem 1.2. (Weak convergence of the lateness function) As R — >• oo, the rescaled 
functions on M 2 defined by Gr((x±, X2)) := ^((L-^iJj L-^ x 2j)) converge to the aug- 
mented Gaussian free field h in the following sense: for each set of test functions 
4>i, . . . , (pk in Hq, the joint law of the inner products (4>j, Gr) converges to the joint 
law of ((f) j, h). 

Our next result addresses the fluctuations from circularity at a fixed time, as 
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Figure 2: Top: Symmetric difference between IDLA cluster Ax(t) at continuous 
time t = 10 5 and the disk of radius \ftpn. Bottom: closeup of a portion of the 
boundary. Sites outside the disk are colored red if they belong to A t m ; sites inside 
the disk are colored blue if they do not belong to Ax{t)- 



illustrated in Figure 2. 

Theorem 1.3. (Fluctuations from circularity) Consider the distribution with point 
masses on M 2 given by 

E t := r _1 ^ (ixeArit) ~ l xeB r )^x/r, (2) 

xez 2 

where r = wt/n. As t —> oo, the Et converge to the restriction of the augmented 
GFF to dB\(0), in the sense that for each set of test functions <j>\, . . . , in H$, 
the joint law of (<pj, E t ) converges to the joint law of &h{4>j, tt) (a Gaussian process 
defined in §1.4)- 

1.4 Main results in general dimensions 

In this section, we will extend Theorem 1.3 to general dimensions and to a range 
of times (instead of a single time). That is, we will try to understand scaling limits 
of the discrepancies of the sort depicted in Figure 2 (interpreted in some sense 
as random distributions) in general dimensions and taken over a range of times. 
However, some caution is in order. By classical results in number theory (see the 
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survey [IKKN04] for their history), the size of B r = B r (0) nZ d is approximately the 
volume of B r (0) — but with errors of order r d ~ 2 (i.e., both 0(r d ~ 2 ) and S7(r d ~ 2 )) 
in all dimensions d > 5. The errors in dimension d = 4 are of order r d ~ 2 times 
logarithmic factors that grow to infinity. It remains a famous open problem in 
number theory to estimate the errors when d £ {2,3}. (When d = 2 this is called 
Gauss's circle problem.) 

These number theoretic results imply that #B r ^) is, as a function of t, much 
more irregular than the size T(t) of the cluster obtained in continuous time internal 
DLA, at least when d > 5. The results also imply that even if points were added 
to At precisely in order of increasing radius, the difference between the radius of At 
and the radius of B r u\(0) would fail to be o(r _1 ) when d > 5 and fail to be 0{r~ l ) 
when d = 4. 

On the other hand, we will see that the kinds of fluctuations that emerge from 
internal DLA randomness are of the order that one would obtain by spreading an 
extra r d l 2 ~ \ft particles over a constant fraction of the spherical boundary, which 
is also what one obtains by changing the radius (along some or all of the boundary) 
by r l ~ d / 2 . This implies that the higher dimensional analog of Theorem 1.3 cannot 
be true exactly the way it is stated if d > 4. Indeed, suppose that we define E t 
analogously to (2) as 

E t = r- d / 2 (Wa t(() - l x en r )S x/r , 

xez d 

and let (ft be a test function that is equal to 1 in a neighborhood of dB\(0). Then 
the results mentioned above imply that 

(E t ,cf>)=r- d / 2 {T(t)-#B r ) 

cannot converge in law to a finite random variable if d > 4. 

It is therefore a challenge to formulate a central limit theorem for the (small) 
fluctuations of internal DLA that is not swamped by these (potentially large) number 
theoretic irregularities. We will see below that this can be achieved by replacing B r 
with different ball approximations (the so-called "divisible sandpiles") that are in 
some sense even "rounder" than the lattice balls themselves. We will also have to 
define and interpret the (augmented) GFF in a particular way. 

Given smooth real- valued functions / and g on M. d , write 

(/,<?) v = / V/(x) • Vg(x)dx. 

Here and below dx denotes Lebesgue measure on M. d . Given a bounded domain 
D C M d , let H(D) be the Hilbert space closure in (-,-)v of the set of smooth 
compactly supported functions on D. We define H = H(M. d ) analogously except 
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that the functions are taken modulo additive constants. The Gaussian free field 
(GFF) is defined formally by 



where the fi are any fixed (•, -)v orthonormal basis for H and the cti are i.i.d. mean 
zero, unit variance normal random variables. (One also defines the GFF on D 
similarly, using H(D) in place of H.) The augmented GFF will be defined similarly 
below, but with a slightly different inner product. 

Since the sum (3) a.s. does not converge within H, one has to think a bit about 
how h is defined. Note that for any fixed f = Y Pifi H, the quantity (h, /)v := 
Yi a ifii /)v = Y a iPi is almost surely finite and has the law of a centered Gaussian 
with variance ||/||v = Y IftP- However, there a.s. exist some functions f E H for 
which the sum does not converge, and (h, -)v cannot be considered as a continuous 
functional on all of H. Rather than try to define (h, /)y for all / £ H, it is often 
more convenient and natural to focus on some subset of / values (with dense span) 
on which / i— > (h, f)y is a.s. a continuous function (in some topology). Here are 
some sample approaches to defining a GFF on D: 

1. h as a random distribution: For each smooth, compactly supported 4>, 
write (h,cj)) := (h, — A _1 0)v, which (by integration by parts) is formally the 
same as J h(x)(p(x)dx. This is almost surely well defined for all such (p and 
makes h a random distribution [Shc07]. (If D = ~R. d and d = 2, one requires 
J <f>(x)dx = 0, so that (h, 4>) is defined independently of the additive constant. 
When d > 2 one may fix the additive constant by requiring that the mean of 
h on B r (0) tends to zero as r — > oo [Shc07].) 

2. h as a random continuous (d + l)-real-parameter function: For each 
e > and x G M d , let h £ (x) denote the mean value of h on dB £ (x). For each 
fixed x, this h E (x) is a Brownian motion in time parameterized by — loge in 
dimension 2, or —e 2 ~ d in higher dimensions [Shc07]. For each fixed e, the 
function h e can be thought of as a regularization of h (a point of view used 
extensively in [DS10]). 

3. h as a family of "distributions" on origin-centered spheres: For each 
polynomial function ip on M. d and each time t, define $/i(V , ) t) to be the integral 
of hip over dB r (0) where B r (0) is the origin-centered ball of volume t. We 
actually lose no generality in requiring if) to be a harmonic polynomial on M. d , 
since the restriction of any polynomial to dB r (0) agrees with the restriction 
of a (unique) harmonic polynomial. 

The last approach turns out to be particularly natural for our purposes. Using 
this approach, we will now give our first definition of the augmented GFF: it is the 
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centered Gaussian function for which 

Cov(* ft (^i,ii),* h (^2,i 2 )) = / ip t (x)ip 2 (x)dx, (4) 

JB r (0) 

for all harmonic polynomials ip\ and ip2, where B r (0) is the ball of volume min{ti, £2}- 
In particular, taking ipi = ip2 = ip, we find that 

Var^hftM)) = / ip{x) 2 dx. (5) 

Though not immediately obvious from the above, we will see in §1.5 that this 
definition is very close to that of the ordinary GFF. Now, for each integer m and 
harmonic polynomial ip, there is a discrete harmonic polynomial V'(m) on (°-e- 
fined precisely in §2.2) that approximates ip in the sense that ip — ipi m \ 1S a P or y- 
nomial of degree at most k — 2, where k is the degree of ip. In particular, if we 
fix ip and limit our attention to x in a fixed bounded subset of then we have 
\ip( m ){x) — ip(x)\ = 0{l/m 2 ). Let Q denote the grid comprised of the edges connect- 
ing nearest neighbor vertices of Z d . (As a set, Q consists of the points in M d with at 
most one non-integer coordinate.) As in [JLSlOa], we extend the definition of ipi m ) 
to Q by linear interpolation. 

Now write 

*2ftM) :=m- d/2 i>(m)(x/m)-tiP {m) (0). (6) 

This random variable measures to what extent the mean value property for the 
discrete harmonic polynomial ip(m) ^ aus f° r ^ ne se ^ ^-T(m d t)- ^ s suc h) it is a way of 
measuring the deviation of ^4;r(m d t) from circularity. 

Theorem 1.4. Let h be the augmented GFF, and $^ as discussed above. Then 
as m —> oo, the random functions converge in law to ^ (w.r.t. the smallest 
topology that makes $ i— > <&(ip,t) continuous for each ip and t). In other words, 
for any finite collection of pairs {ip\,t\), . . . , {ipk,tk)> the joint law of the <&™(ipi,ti) 
converges to the joint law of the <&h(ipi,ti). 

Remark. Theorem 1.4 does not really address the discrepancies between B r and At 
(which, as we remarked, can be very large, in particular in the case that ip is a 
constant function). However, it can be interpreted as a measure of the discrepancy 
between At and the so-called divisible sandpile, which is a function wt : Z rf — )• [0, 1] 
defined for all t > 0. The quantity wt(x) represents the amount of mass that 
ends up at x if one begins with t units of mass at the origin and then "spreads" 
the mass around according to certain rules that ensure that the final amount of 
mass at each site is at most one. For fixed x, the quantity wt(x) is a continuously 
increasing function of t, and moreover there exists a constant c depending only on 
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the dimension d, such that Wt(x) = 1 if |x| < r(t) — c and Wt(x) = if \x\ > r(t) + c 
[LP09]. An important property of Wt(x) is that for any discrete harmonic function 
ip on Z d we have Yl x ^z d w t(x)tp(x) = 0. It is natural to replace (2) with 

Ek := r~ d / 2 £ (W T(t) ~ M^x/r, (7) 



and interpret Theorem 1.4 as a statement about the distributional limit of Et- 

Remark. Even with the replacement above, Theorem 1.4 differs from Theorem 1.3, 
since it requires that we use only harmonic polynomial test functions tp and also 
requires that we replace them with approximations ^( m ) on the discrete level. It 
is natural to ask, in general dimensions, what happens when we try to modify the 
statement of Theorem 1.4 (interpreted as a sort of distributional limit statement 
for (7)) to make it read like the distributional convergence statement of Theorem 
1.3. We will discuss this in more detail in §3.4, but we can summarize the situation 
roughly as follows: 



Modification 


When it matters 


Replacing wt in (7) with le r 


No effect when d = 2. 
Invalidates result when d > 3. 


Keeping wt in (7) but 
Replacing ipr m ) with ip 


No effect if d € {2,3,4,5}. 
Unclear if d > 5. 


Keeping wt in (7) but 
Replacing tpf m ) with general 
smooth test function 6. 


No effect if d € {2,3}. 

Probably invalidates result if d > 3. 



The restriction to harmonic (as opposed to a more general test function <j>) 
seems to be necessary in large dimensions because otherwise the derivative of the 
test function along dB\(0) appears to have a non-trivial effect on (7) (see §3.4). This 
is because (7) has a lot of positive mass just outside of the unit sphere and a lot of 
negative mass just inside the unit sphere. It may be possible to formulate a version of 
Theorem 1.4 (involving some modification of the "mean shape" described by wt) that 
uses test functions that depend only on 9 in a neighborhood of the sphere (instead of 
using only harmonic test functions), but we will not address this point here. Deciding 
whether Theorem 1.2 as stated extends to higher dimensions requires some number 
theoretic understanding of the extent to which the discrepancies between wt and 
(as well as the errors that come from replacing a ^(m) with a smooth test function 
(f>) average out when one integrates over a range of times. We will not address these 
points here either. 
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1.5 Comparing the GFF and the augmented GFF 

We may write a general vector in W 1 as rO where r £ [0,oo) and 9 £ S^ 1 := &B X (0). 
We write the Laplacian in spherical coordinates as 

A = r ^dd d ^d_ 2 

dr dr b w 

A polynomial / G M[a;i, . . . , x^] is called harmonic if A/ is the zero polynomial. 
Let Vi denote the space of all homogenous harmonic polynomials in . . . , Xd] of 
degree £, and let Hp denote the space of functions on S^ 1 obtained by restriction 
from Ve. If / E He, then we can write f(rO) = g(6)r e for a function g £ i^, and 
setting (8) to zero at r = 1 yields 

A 5d -isr = -£(£ + d-2)g, 

i.e., (7 is an eigenfunction of Agn-i with eigenvalue —£{£ + d — 2). Note that (8) 
continues to be zero if we replace £ with the negative number £' := — {d — 2) — £, 
since the expression — £(£ + d — 2) is unchanged by replacing £ with £'. Thus, g(9)r 
is also harmonic on M. d \ {0}. 

Now, suppose that g is normalized so that 

/ g{0fd6 = 1. (9) 

By scaling, the integral of f 2 over <9-Br(0) is thus given by R d ~ 1 R 2e . The L 2 norm 
on all of -Br(O) is then given by 



Br(0) 



■R T>d+2i 

J{zfdz = I r*-W = J^. (10) 



A standard identity states that the Dirichlet energy of g, as a function on S ^, 
is given by the L 2 inner product (— Ag, g) = £(£+d— 2). The square of ||V/|| is given 
by the square of its component along S^ -1 plus the square of its radial component. 
We thus find that the Dirichlet energy of / on -Br(O) is given by 



r-R f-R 

\\Vf{z)\\ 2 dz = £{£ + d-2) / r d-i r 2(e-i) dr + / r d-i r 2(i-i)^2 dr 



b r (0) Jo 



£{£ + d — 2) 2 £+d-2 , ^ 2 p2l+d-2 

2£ + d-2 2£ + d-2 

2£ 2 + (d-2)£ 2£+d _ 2 
2£+{d-2) 
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Now suppose that we fix the value of / on dBfi(0) as above but harmonically 
extend it outside of Br(0) by writing f(r0) = R l ~ l g(0)r e for r > R. Then the 
Dirichlet energy of / outside of -Br(O) can be computed as 

poo poo 

RW-V)l(t + d-2) / r^rW-Vdr + RW-^ / r^r^'^ie'fdr, 
Jr Jr 

which simplifies to 

e + i(d - 2) + (£T R 2i +d -2 = p + gj - 2) + (i + (d - 2)) 2 R2e+d _ 2 



2£' + (d-2) 2(-£- (d-2)) + (d-2) 

-i? 



2£ 2 + 3£(d-2) + (d-2) 2 o2 , +d _ 



-21 - (d - 2) 
= (^ + d_2) J R 2 ^- 2 . 

Combining the inside and outside computations in the case R = 1, we find that 
the harmonic extension / of the function given by g on S^ -1 has Dirichlet energy 
2£ + (d — 2). If we decompose the GFF into an orthonormal basis that includes this 
/, we find that the component of / is a centered Gaussian with variance 2 £+(d-2) • 
If we replace / with the harmonic extension of g(R~ 1 6) (defined on <9.Br(0)), then 
by scaling the corresponding variance becomes 2 e+(d-2) R 2 d - 

Now in the augmented GFF the variance is instead given by (10), which amounts 
to replacing 2 e+(d-2) w ^ ri w+d' Considering the component of g(R^ 1 9) in a basis 
expansion the space of functions on dBji(0) requires us to divide (10) by R 2£ (to 
account for the scaling of /) and by (R^ 1 ) 2 (to account for the larger integration 
area), so that we again obtain a variance of 2 jj^R 2 ^ d for the augmented GFF, 

versus 2l+ l d ^ 2) R 2 ~ d for the GFF. 

In some ways, the augmented GFF is very similar to the ordinary GFF: when 
we restrict attention to an origin-centered annulus, it is possible to construct inde- 
pendent Gaussian random distributions h\, h.2, and /13 such that h\ has the law of 
a constant multiple of the GFF, h\ + hi has the law of the augmented GFF, and 
hi + hi + /13 has the law of the ordinary GFF. 

In light of Theorem 1.3, the following implies that (up to absolute continuity) 
the scaling limit of fixed-time At fluctuations can be described by the GFF itself. 

Proposition 1.5. When d = 2, the law v of the restriction of the GFF to the unit 
circle (modulo additive constant) is absolutely continuous w.r.t. the law fi of the 
restriction of the augmented GFF restricted to the unit circle. 

Proof. The relative entropy of a Gaussian of density e~ x2 ^ 2 with respect to a Gaus- 
sian of density a~ 1 e~ x l^ 2a is given by 

F(a) = J e~ x2 / 2 ({a~ 2 - l)x 2 /2 + log a) dx = (a' 2 - l)/2 + log a. 
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Note that F'(a) = — <r -3 + a^ 1 , and in particular -F'(l) = 0. Thus the relative 
entropy of a centered Gaussian of variance 1 with respect to a centered Gaussian 
of variance 1 + a is 0(a 2 ). This implies that the relative entropy of \x with respect 
to v — restricted to the jth component ctj — is 0(j~ 2 ). The same holds for the 
relative entropy of v with respect to \i. Because the otj are independent in both /i 
and the relative entropy of one of \x and v with respect to the other is the sum 
of the relative entropies of the individual components, and this sum is finite. □ 

2 General dimension 

2.1 FKG inequality: Proof of Theorem 1.1 

We recall that increasing functions of a Poisson point process are non-negatively 
correlated [GK97]. (This is easily derived from the more well known statement 
| 171] that increasing functions of independent Bernoulli random variables are 
non-negatively correlated.) Let /i be the simple random walk probability measure 
on the space Q' of walks W beginning at the origin. Then the randomness for 
internal DLA is given by a rate-one Poisson point process on fi x v where v is 
Lebesgue measure on [0, oo). A realization of this process is a random collection of 
points in Q' x [0, oo). It is easy to see (for example, using the abelian property of 
internal DLA discovered by Diaconis and Fulton [DF91]) that adding an additional 
point (w, s) increases the value of A T u\ for all times t. The A T n\ are hence increasing 
functions of the Poisson point process, and are non-negatively correlated. Since F 
and G are increasing functions of the A T u\ , they are also increasing functions of the 
point process — and are thus non-negatively correlated. 

2.2 Discrete harmonic polynomials 

Let i/)(xi, . . . , Xd) be a polynomial that is harmonic on that is 

In this section we give a recipe for constructing a polynomial ipi that closely ap- 
proximates tp and is discrete harmonic on Z d , that is, 

d 

1=1 

where 

£> 4 2 Vi = i/)i(x + e f ) - 2^i (x) +ij>i(x- ej) 

is the symmetric second difference in direction ej. The construction described below 
is nearly the same as the one given by Lovasz in [Lov04] , except that we have tweaked 
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it in order to obtain a smaller error term: if ip has degree k, then ip — ip\ has degree 
at most k — 2 instead of k — 1. Discrete harmonic polynomials have been studied 
classically, primarily in two variables: see for example Duffin [Duf56], who gives a 
construction based on discrete contour integration. 
Consider the linear map 

E : R[x 1 , ...,x d }^ R[xi, ...,x d ] 

defined on monomials by 

Six* 1 ■ ■ ■ x k /) = P kl ( Xl ) ■ ■ ■ P kd (x d ) 

where we define 

(fc-l)/2 

Pk(x)= n ( x+ ti- 

i=-(fc-l)/2 

Lemma 2.1. If ip E M.[x\, . . . , x d ] is a polynomial of degree k that is harmonic on 
M. d , then the polynomial ipi = H(^) is discrete harmonic on Z rf , and ip — tpi is a 
polynomial of degree at most k — 2. 

Proof. An easy calculation shows that 

D 2 P k = k(k - l)P fc _ 2 

from which we see that 

If tp is harmonic, then the right side vanishes when summed over i = 1, . . . , d, which 
shows that E[ip] is discrete harmonic. 

Note that P k {x) is even for k even and odd for k odd. In particular, P k {x) — x k 
has degree at most k — 2, which implies that tp — ipi has degree at most k — 2. □ 

To obtain a discrete harmonic polynomial ipt m ) on the lattice ^^ d , we set 

fp(m){ x ) '■= rrr k ipi{mx), 

where k is the degree of tp. 

2.3 General-dimensional CLT: Proof of Theorem 1.4 

Proof of Theorem 1.4. For each fixed tp with V'(O) = the process M(t) := 
&™(ip,t) is a martingale in t. Each time a new particle is added, we can imagine 
that it performs Brownian motion on the grid (instead of a simple random walk), 
which turns M into a continuous martingale, as in [JLSlOa]. By the martingale 
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representation theorem (see [RY05, Theorem V.1.6]), we can write M(t) = /3(s m (t)), 
where (3 is a standard Brownian motion and 



n-l 

s m (t):= limsup V(M(t m ) - M(U)) 2 
0=t <ti<-<t„=t . „ 
*»|— >0 

is the quadratic variation of M on the interval [0, t]. To show that ^^{ip, t) converges 
in law as m — > oo to a Gaussian with variance V := J B , Q \ tp(x) 2 dx, it suffices to 

show that for fixed t the random variable s m (t) converges in law to V. 

By standard Riemann integration and the At fluctuation bounds in [JLSlOa, 
JLSlOb] (the weaker bounds of [LBG92] would also suffice here) we know that 
m~ d J2xeA d ^{m)^ l m ) 2 converges in law to f B ip(x) 2 dx as m — > oo. Thus it 
suffices to show that 

s m (t)-m- d V i{; {m) {x/m) 2 (11) 



Km" 



converges in law to zero. This expression is actually a martingale in t. Its expected 
square is the sum of the expectations of the squares of its increments, each of which 
is 0(m~ 2d ). The overall expected square of (11) is thus 0(m~~ d ), which indeed tends 
to zero. 

Recall (6) and note that if we replace t with T(t), this does not change the 
convergence in law of s m (t) when ^(0) = 0. However, when V'(O) 7^ 0, it introduces 
an asymptotically independent source or randomness which scales to a Gaussian of 
variance ip(0) 2 t (simply by the central limit theorem for the Poisson point process), 
and hence (5) remains correct in this case. 

Similarly, suppose we are given = to < t\ < ti < . . . < t£ and functions 
■01, ipz, . . . tpi. The same argument as above, using the martingale in t, 



m 



-d/2 



^2 ^j,(m)(s/ra)-%,(m)(0)> 



j 1 Ar(m d (tAt,-)) 



implies that X^i=i C0j> tj) converges in law to a Gaussian with variance 



E 



/ y^fpifc) dx. 

JB„n.\\B„u. .1 \ / 



j = l jB r(t j )\B r(tj _ 1 ) \ i= j 

The theorem now follows from a standard fact about Gaussian random variables on 
a finite dimensional vector spaces (proved using characteristic functions): namely, a 
sequence of random variables on a vector space converges in law to a multivariate 
Gaussian if and only if all of the one-dimensional projections converge. The law of h 
is determined by the fact that it is a centered Gaussian with covariance given by 
(4). " □ 
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3 Dimension two 

3.1 Two dimensional CLT: Proof of Theorem 1.2 

Recall that At for t G Z + denotes the discrete-time IDLA cluster with exactly t sites, 
and At = A T (t) fo r t E M+ denotes the continuous-time cluster whose cardinality is 
Poisson-distributed with mean t. 
Define 

F (t) := inf{t :z€A t } 

and 

L (z) := y/F Q (z)/n-\z\. 
Fix N < oo, and consider a test function of the form 

\k\<N 

where the a k are smooth functions supported in an interval < tq < r < r\ < oo. 
We will assume, furthermore, that <p is real-valued. That is, the complex numbers 
a k satisfy 

a_fc(r) = a k {r) 

Theorem 3.1. As R — > oo, 

& E L (Rz)^l ^ N(0,V ) 



in law, where 



/*oo /*oo 



0<|fc|<7V 



dr 
r 



2 dp 



It follows from Lemma 3.2 below (with q = \k\ + 1, y = logr), that Theorem 3.1 
can be interpreted as saying that Lq(Rz) tends weakly to the Gaussian random 
distribution associated to the Hilbert space H\ r with norm 

. /-oo J 



l»7llo 

0<|fc|<oo 



V 2tt / [\rd r r, k \ 2 + {\k\ + lf\r, k \ 



where 



1 



2- 



o 



and 770 (r) = 0. (The subscript nr means nonradial: H\ r is the orthogonal comple- 
ment of radial functions in the Sobolev space H 1 .) 
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Lemma 3.2. Let q > and let ip be a real valued function on R. Denote 

/oo 
i>(y)f{y)dy 
-oo 

where the supremum is over real-valued f, compactly supported and subject to the 
constraint 



/oo 
(f'(y) 2 + q 2 f(y) 2 )dy<i. 
-oo 



Then 



i/j(y)e q ( s ~ y) dy 
Proof. In the case q = 0, replace / in / ?/>/ dy with 



ds. 



/(y) 



change order of integration and apply the Cauchy-Schwarz inequality. For the case 
q > 0, multiply by the appropriate factors e qy and e 2qy to deduce this from the case 
q = 0. □ 

If we use Ay and corresponding functions F(z) and L(z), then the ao coefficient 
figures in the limit formula as follows. 



Theorem 3.3. As R — > oo, 



i? 2 



7V(0,F) 



ze(z+jZ)/_R 



in Zaw, where 



/OO /'OO 
/ a fe (r)(p/r)l fc l +1 



dr 
r 



dp 



Theorem 3.3 is a restatement of Theorem 1.2. (As in the proof of Theorem 1.4, 
the convergence in law of all one-dimensional projections to the appropriate normal 
random variables implies the corresponding result for the joint distribution of any 
finite collection of such projections.) Theorem 3.3 says that L(Rz) tends weakly to 
a Gaussian distribution for the Hilbert space H 1 with the norm 



°° /-OO 

V 2tt / [|r^»/ fc | 2 + (|A;| + l) 2 |^| 2 ] 
Jo 



k=—oo 



dr 
r 
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By way of comparison, the usual Gaussian free field is the one associated to the 
Dirichlet norm 



/* /*00 J 

/ \V V \ 2 dxdy = V 2vr / [\rd rm \ 2 + k 2 \ m \ 2 f r 
Jr 2 ' Jo 



r 



k=—oo 

Comparing these two norms, we see that the second term in \\rj\\ 2 has an additional 
+1, hence our choice of the term "augmented Gaussian free field." As derived in 
§1.5, this +1 results in a smaller variance 2E+d^ d * n eacn spherical mode of degree 
£ of the augmented GFF, as compared to 2 £+d-2 ^ 2 ~ d ^ or ^ e usual GFF. The surface 
area of the sphere is implicit in the normalization (9), and is accounted for here in 
the factors 2ir above. 



To prove Theorem 3.1, write 

Lo(z) = -?= f°(l - U t )* 1/2 y - f°(l - U z? <t)t 

W n Jo t 2 V 7T Jo 

/ (K\z\ 2 <t - i-Ajt '"— 
Jo 



2Vvr J " - l " ty t 
Let po(z) = 1, and for k > 1 let Pk(z) = qk(z) — 9fc(0), where 

q k (z) =E[z k ] 

is the discrete harmonic polynomial associated to z k = {x + iy) k as described in 
§2.2. The sequence pk begins 

i 2 3—4 - 

.L ^ Z j Z •) Z ~~ Z , Z Z Z j ... ■ 

For instance, to compute ps, we expand 

3 _ „,3 o_,2 , • [o„2„, 31 



Z = X 



3xy 2 + i [3x 2 y — y 3 



and apply E to each monomial, obtaining 
Pz{z) = (x-l)x(x+l)-3x(y-~)(y + ~) + i 



3(x--)(x + -)y-(y-l)y(y + l) 



which simplifies to z 3 — \z. One readily checks that this defines a discrete harmonic 
function on Z + iL. (In fact, z 3 is itself discrete harmonic, but z k is not for k > 4.) 
To define pk for negative k, we set p-k(z) = Pk(z)- 
Define 

N 

i>(z,t,R)= J2 MVtF^)Pk(z)(VtFr lkl 

k=-N 

and 

Mz, t, R) = ip{z, t, R) - ao{Vt/KR 2 ) 
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Lemma 3.4. If c\R 2 <t< C2R 2 and \ \z\ — y/t/iv\ < C\ogR, then 

\^(z,t,R)-4>{z/R)\<C{\ogR)/R 

This lemma follows easily from the fact that the coefficients are smooth and 
the bound \pk(z) - z k \ < C|^|l fe l -1 . 

3.2 Van der Corput bounds 
Lemma 3.5. (Van der Corput) 

(a) \#{z £ Z + iZ : ir\z\ 2 < t} - t\ < Ct 1 / 3 . 

(b) For k>l, 



-k/2 



^ zk ^n\z\ 2 <t 

zeZ+iZ 



< Ct 1 ' 3 . 



(c) For k > 1, 



,-k/2 



\z\ 2 <t 



< Ct l '\ 



Part (a) of this lemma was proved by van der Corput in the 1920s (See [GS10], 
Theorem 87 p. 484). Part (b) follows from the same method, as proved below. Part 
(c) follows from part (b) and the stronger estimate of Lemma 2.1, \pk(z) — z k \ < 
C\z\ k ~ 2 for k > 2 (and Pl {z) -z = 0). 

We prove part (b) in all dimensions. Let Pk be a harmonic polynomial on M rf of 
homogeneous of degree k. Normalize so that 

max |Pfc(x)| = 1 

x£B 

where B is the unit ball. In this discussion k will be fixed and the constants are 
allowed to depend on k. 

We are going to show that for k > 1, 



\x\<R,xGZ d 



< R 



-l-a 



where 



a = 1 - 2/(d+ 1) 

For d = 2, a = 1/3, and R d R- 1 ~ a = R 2 / 3 ps t 1 / 3 . This is the claim of part (b). 



The van der Corput theorem is the case k = 0. It says 
(l/R d ) |x| < r\ - vol (|x| < R) 



< R 



-l-a 
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Let e = l/R a . 

Consider p a smooth, radial function on R d with integral 1 supported in the unit 
ball. Then define x = Ifi characteristic function of the unit ball. Denote 

p e (x) = e' d p(x/e), XR (x) = R- d X {x/R) 

Then 

J2(XR*Pe(x)-XR(x))Pk(x)/R< 

xez d 

This is because \R * Pt( x ) ~ Xr( x ) is nonzero only in the annulus of width 2e around 
\x\ = R in which (by the van der Corput bound) there are 0(R d ~ 1 e) lattice points. 
The Poisson summation formula implies 

£ xr * Pe (x)P k (x)/R k = Y, \xr(OMQ] * h(0/R k 

in the sense of distributions. The Fourier transform of a polynomial is a derivative 
of the delta function, Pk(0 = Pk(id^)S(^). Because k > 1 and Pk(x) is harmonic, 
its average with repect to any radial function is zero. This is expressed in the dual 
variable as the fact that when £ = 0, 

Pk(ids)(x R (OPe(0])=0 

So we our sum equals 

E lxR(Ope(0]*Pk(0/R k 

Next look at 

xr(0 = x(R0 

P k (idt) X (R0 =R k I Pk{x)e- mx <dx 

J\x\<\ 

All the terms in which fewer derivatives fall on xr an d more fall on p e give much 
smaller expressions: the factor R corresponding to each such differentiation is re- 
placed by an e. 

The asymptotics of this oscillatory integral above are well known. For any fixed 
polynomial P they are of the same order of magnitude as for P = 1, namely 

\P k (id0x(m\/R k <C k \R£\- (d+1)/2 

This is proved by the method of stationary phase and can also be derived from well 
known asymptotics of Bessel functions. 
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It follows that our sum is majorized by (replacing the letter d by n so that it 
does not get mixed up with the differential dr) 



s: 



(Rr) 



>+!)_ 



„rc-l 



dr 



(1 + er) 



N 



1/6 



{Rr 



-(n+l) r n 



dr 



1 r 

jj-(n+l)/2 e -(n-l)/2 



3.3 The main part of the proof 

Denote 



L E L o(^)- 



R 2 



z£(Z+iZ)/R 



Applying the formula above for Lq, 

x R= £ 



\z\ 2 <t 



1 



E (i 



t|z| 2 <£ 



H z / R ) 4 l/2 dt 

\z\ 2 t 

4>(z,t,R) 1/2 dt 

— — 1 — + h R . 

t/TT t 



To estimate the error term E R , note first that the coefficients a& are supported 
in a fixed annulus, the integrand above is supported in the range c\R 2 < t < 
C2R 2 . Furthermore, by [JLSlOa], there is an absolute constant C such that for all 
sufficiently large R and all t in this range, the difference l % \ z \2 <t — \& t is supported 
on the set of z £ Z 2 such that \\z\ — \/t/ir\ < ClogR. Thus 

\l v \ z \*< t -lA t \<KR]ogR. 



Moreover, Lemma 3.4 applies and 

-c 2 R 2 

' C1 R 2 

Next, Lemma 3.5(a) says (since #A t = t) 



\E R \ < C / ORlogi^^V 1 / 2 - = 0((logR)/R). 

J a R 2 K t 



£ 1 



ir\z\ 2 <t 



< ct 1 / 3 . 
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Thus replacing ip by -00 gives an additional error of size at most 

rc 2 R 2 Ji 

C / t Vh~ l ' 2 - = 0{R- 1 ^). 

In all, 

/•7r r°° fit 



12) 



For s = 0, 1, . . . , consider the process 



M(s) 



i/ 2 dt 



Note that M(s) —> Xr as s — > oo. Note also that Lemma 3.5(c) implies 

M(0) = 0(iT 1/3 ). 

Because are discrete harmonic and Pfc(0) = for all k ^ 0, M(s) — M(0) is 
a martingale. It remains to show that M(s) — M(0) — > N(0, Vq) in law. As 
outlined below, this will follow from the martingale central limit theorem (see, e.g., 
[Bro71, HH80] or [Dur95, p. 414]). 

For sufficiently large R, the difference M{s + 1) — M(s) is nonzero only for s in 
the range c\R 2 < s < c 2 R 2 \ and \Fq(z) — tt\z\ 2 \ < CRlogR. We now show that this 
implies 

\M(s + l) - M(s)\ 2 = 0(1/R 2 ) (13) 

and 

00 

\M(s + 1) - M(s)\ 2 = V + 0((logR)/R) (14) 

s=0 

so that the martingale central limit theorem applies. 
To prove (13), observe that 



/tt f°° 

M( s + 1) - M(s) = ~ M*, t> R )t 

1 JFn(z) 



Fo(z) 1 

where z is the (s+l)th point of A t . Then \z\ < ^/T/tt+K log R implies \pk(z)\(t/-ir)-\ k ^ 2 < 
C, and hence 

\Mz,t,R)\ < C 
Recalling that ipg = unless c\R 2 < t < C2R 2 , we have 



[■C2R 2 f if 
\M{s + l)-M{s)\<C r l/2 - =0(1/R) 

J a B? t 
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which confirms (13). 

Because A t fills the lattice Z + iL as t — > oo, we have 



Y^\M(s + l)- M(s) 



s=0 



>7T 



Fo(z) 



a k {^i/^)p k {z){t/^)-\ k \l 2 r 1 



/2 



0<|fc|<JV 



We prove (14) in three steps: replace p k (z) by z fc (or z^ if A; < 0); replace the lower 
limit Fq(z) by vr|z| 2 ; replace the sum of z over lattice sites with the integral with 
respect to Lebesgue measure in the complex z-plane. 

We begin the proof of (14) by noting that the error term introduced by replacing 
p k with z k is 

\p k {z) - z k \(t/7r)-W < C^- 1 = 0(1/R 2 ) 



In the integral this is majorized by 

C2R \-i/2dt 
ciH 2 * 



C2-R 2 1 rff 



Since there are 0(R 2 ) such terms, this change contributes order R 2 /R = 1/R 2 to 
the sum. 

Next, we change the lower limit from Fq(z) to vr | 2 . Since |i 7 b( 2; ) — ^M 2 ! < 
CRlogR, the integral inside | • • • | 2 is changed by 



0{{\ogR)/R 2 



"n\z\ 2 j + 

l Cl R?<c 2 R 2t 
'Fo(z) * 

Thus the change in the whole expression is majorized by the order of the cross term 

(1/ R)(log R) / R 2 = (log R)/R 3 

Again there are R 2 terms in the sum over z, so the sum of the errors is 0((log R)/R). 
Lastly, we replace the value at each site zq by the integral 



a k {^/t/TTTP)r k e ik9 {t/TT)-WI 2 r l 



/2 



dt 



0<|fc|<AT 



rdrdO 



where Q ZQ is the unit square centered at zq and z = re . Because the square has 
area 1, the term in the lattice sum is the same as this integral with z = re %e replaced 
by zq at each occurrence. Since \z — zq\ < v2, 



\z k -4\ < 4fc (M 



\fc-l 
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After we divide by (y / t/7r) k , the order of error is 1/R. Adding all the errors con- 
tributes at most order 1/R to the sum. We must also take into account the change in 
the lower limit of the integral, 7r|zo| 2 is replaced by tt\z\ 2 = irr 2 . Since \z — Zq\ < y/2, 



< V2(\z\ + \zo\) < CR 



Recall that in the previous step we previously changed the lower limit by 0(R log R). 
Thus by the same argument, this smaller change gives rise to an error of order 1/R 
in the sum over zq. 

The proof of (14) is now reduced to evaluating 



2-7T pOO 

J 



rdrdO 



0<\k\<N 

Integrating in 9 and changing variables from r to p = r /R, 

V / / a k {^J^){Rp)\ k \ + \t/^)-\ k \lh^\ 

lnp 2 R 2 * 



0<|fcl<iV 







dp 
P 



Then change variables from t to to r = \ft/irR 2 to obtain 



2n £ 

0<|fc|<Af 



a k (r)(p/rp +1 



dr 



dp 
P 



V . 



This ends the proof of Theorem 3.1. 

The proof of Theorem 3.3 follows the same idea. We replace At by the Poisson 
time region At (for T = T(t)), and we need to find the limit as R — > oo of 



V/ (t-#A T )a (^7^r 1/2 

1 Jo 



dt 



poo 

+ ^ ^N 2 ^* ~ 1 A T )ipo(z,t,R)r 



1/2 



The error terms in the estimation showing this quantity is within 0(R x / 3 ) of 



R 2 



are nearly the same as in the previous proof. We describe briefly the differences. 
The difference between Poisson time and ordinary counting is 



\#A T - #A t \ = \#A T -t\< Ct 1 ' 2 log* = O(RlogR) almost surely 
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if t ~ R 2 . It follows that for \z\ « i?, 

|.F(z) — vr|z| 2 | = O(RlogR) almost surely 

as in the previous proof for i*b(z)- Further errors are also controlled since we then 
have the estimate analogous to the one above for At, namely 



E I l K\z\*<t ~ U T I < CR log R 



We consider the continuous time martingale 



M(s) 



(s A t - #A T{sAt) )ao(y/t/nR 2 )t 



-1/2 dt 



dt 
T 



Instead of using the martingale central limit theorem, we use the martingale repre- 
sentation theorem. This says that the martingale M(s) when reparameterized by 
its quadratic variation has the same law as Brownian motion. We must show that 
almost surely the quadratic variation of M on < s < 00 is 

lim E ((M(s + e) - M{s)) 2 \A T{s) ) / e 



1 

2^ 



2tt 



-1/2 



a k (VTF&y kd )(s/tf\/ 2 t-' 

\k\<N 



dt 



dO 



+ 0{R- 1 / 3 ) 



Integrating with respect to s gives the quadratic variation V + 0(R 1 ^ 3 ) after a 
suitable change of variable as in the previous proof. 



3.4 Fixed time fluctuations: Proof of Theorem 1.3 

Theorem 1.3 follows almost immediately from the d = 2 case of Theorem 1.4 and 
the estimates above. Consider (0, Et) where E% is as in (7). What happens if we 
replace (f> with a function that is discrete harmonic on the rescaled mesh mT l 7h d 
within a log m/m neighborhood of i?i(0)? Clearly, if <j> is smooth, we will have 
(/) — (/) = 0(?n _1 log m). Since there are at most 0{m d ~ l log m) non-zero terms in 
(7), the discrepancy in 

(<t>,E t ) - (4>,E t ) = 0(m~ d / 2 m d '' 1 (mr 1 logm)logm) = 0(m d/2 ~ 2 (log m) 2 ), (15) 

which tends to zero as long as d £ {2,3}. 
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The fact that replacing Et with Et has a negligible effect follows from the above 
estimates when d = 2. This may also hold when d = 3, but we will not prove it 
here. Instead we remark that Theorem 1.3 holds in three dimensions provided that 
we replace (2) with (7), and that the theorem as stated probably fails in higher 
dimensions even if we make a such a replacement. The reason is that (7) is positive 
at points slightly outside of B r (or outside of the support of wt) and negative at 
points slightly inside. If we replace a discrete harmonic polynomial t/j with a function 
that agrees with ip on -Bi(O) but has a different derivative along portions of dB\(0), 
this may produce a non-trivial effect (by the discussion above) when d > 4. 

Finally, we note that replacing ip m by ip introduces an error of order m~ 2 , and 
the same argument as above gives 

(ip,E t ) - (4> m ,Et) = 0(m~ d/2 m d - 1 m- 2 logm) = 0(m d/2 ~ 3 (log m)) , (16) 

which tends to zero when d G {2, 3, 4, 5}. 
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